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Abstract 

It is proved that for general, not necessarily periodic quasi one dimensional systems, the 
band position operator corresponding to an isolated part of the energy spectrum has discrete 
spectrum and its eigenfunctions have the same spatial localization as the corresponding spec- 
tral projection. As a consequence, an eigenbasis of the band position operator provides a basis 
of optimally localized (generalized) Wannier functions for quasi one dimensional systems, and 
this proves the strong Marzari-Vanderbilt conjecture. If the system has some translation sym- 
metries (e.g. usual translations, screw transformations), they are "inherited" by the Wannier 
basis. 

1 Introduction 

Wannier functions (WF) were introduced by Wannier in 1937 [1] as bases in subspaces of states 
corresponding to energy bands in solids, bases consisting of exponentially localized functions (local- 
ized orbitals). For periodic crystals they are defined as Fourier transform of Bloch functions of the 
corresponding bands. Since then WF proved to be a key tool in quantum theory of solids as they 
provide a tight binding description of the electronic band structure of solids. At the conceptual 
level they lay at the foundation of all effective mass type theories e.g the famous Peierls-Onsager 
substitution describing the dynamics of Bloch electrons in the presence of an external magnetic 
field (see e.g.[5]and references therein). At the quantitative level, especially after the seminal pa- 
per by Marzari and Vanderbilt 3 , WF become an effective tool in ab initio computational studies 
of electronic properties of materials. Moreover during the last decades WF proved to be an es- 
sential ingredient in the study of low dimensional nanostructures such as linear chains of atoms, 
nanowires, nanotubes etc (see e.g. [4], [5]). In particular WF are essential for most formulations 
of transport phenomena using real space Green's function method based on Landauer-Biittiker 
formalism both at rigorous [B] and computational levels [7],|I]. 

A few remarks are in order here. The first one is that realistic low dimensional systems are 
not strictly one (two) dimensional but rather quasi one (two) dimensional and one has to take 
into account the (restricted) motion along perpendicular directions. This adds specific features as 
for example the screw symmetry in nanotubes and nanowires absent in strictly one dimensional 
systems. The second one is that realistic systems, due to the presence of defects, boundaries, 
randomness etc, do not have usually full translation symmetry and this ask for a theory of WF 
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not based on Bloch formalism. Finally let us remind that contrary to a widespread opinion (see 
e.g. the discussion in j^j) that WF always exist for isolated band in solids this is not true. More 
precisely, in more than one dimension there are subtle topological obstructions and these are related 
to the QHE [5], [S], [TO]: a band for which WF are known to exist gives no contribution to the 
quantum Hall current. It is then crucial to have rigorous proofs of the existence of exponentially 
localized WF. 

For one dimensional periodic systems the existence of exponentially localized WF has been 
proved by Kohn in his classic paper [llj about analytic structure of Bloch functions. An extension 
of Kohn analysis to quasi one dimensional systems has been done recently by Prodan |12j . As for 
higher dimensions it was known since the work by des Cloizeaux [13] [14j that there are obstructions 
to the existence of exponentially localized WF and that these obstructions are of topological origin 
(more precisely as explicitly stated in [TS] these obstructions are connected to the topology of a 
vector bundle of orthogonal projections). The fact that for simple bands of time reversal invariant 
systems the obstructions are absent was proved by des Cloizeaux [T5] [T3] under the additional 
condition of the existence of centre of inversion and by Nenciu [15] in the general case. While 
the proofs in [13] [14], [15] did not use the vector bundle theory it was suggested in [2] . [16] that 
the characteristic classes theory in combination with some deep results in the theory of analytic 
functions of several complex variables (Oka principle) can be used to give alternative proof of 
the above results and to extend them to composite bands of time reversal symmetric systems. 
This has been substantiated recently in [T7], [T0J where the existence of exponentially localized 
Wannicr functions has been proved for composite bands of time reversal symmetric systems in 
two and three dimensions settling in the affirmative a long standing conjecture. In conclusion the 
situation is satisfactory as far as periodic time reversal symmetric Hamiltonians are considered (as 
already mentioned for Hamiltonians which are not time reversal symmetric exponentially localized 
Wannier functions might not exists). 

As already said above both the theory and applications of Wannier functions boosted since 
Marzari and Vanderbilt [3], introduced studied and proposed methods to compute the so called 
maximally localized Wannier functions (MLWFs) defined by the fact that they minimize the 
position mean square deviation. It was conjectured in [3] that they can be chosen to be real 
functions and that they have "optimal" exponential localization in the sense that they have the 
same exponential localization as the integral kernel of the projection operator of the corresponding 
band. MLWFs proved to be an invaluable tool in the theory of electronic properties of periodic 
media especially in the modern theory of electronic polarizability (see e.g. [18j and references 
therein). 

In the one dimensional case the theory of MLWFs is much more developed. It is known [3] 
that MLWFs are identical to the eigenfunctions of the " band position" operator and then they are 
unique (up to uninteresting phases) and can be chosen to be real functions. Moreover the phases 
of the corresponding Bloch functions are related to the parallel transport procedure [3], [T9] . 
Recently a detailed study of Wannier functions, including their exponential decay, emphasizing 
the difference between the cases with and without inversion symmetry appeared in |20j . In the 
same paper there are pointed out situations in which the Wannicr functions could decay slower 
than the kernel of the projector, which shows that choosing the optimal phase is not a trivial task. 
Our results show that that by choosing the right phase one must always obtain an optimal decay. 

Motivated by the great interest in nonperiodic structures much effort has been devoted to 
extend the results about existence of exponentially localized bases for isolated bands in nonperiodic 
systems. The basic difficulty stems from the fact that for nonperiodic systems one cannot define 
Wannier functions as Fourier transforms of the Bloch functions. One way out of the difficulty 
is to start from the periodic case or tight-binding limit where the Wannier functions are known 
to exist and and use perturbation or "continuity" arguments. The basic idea is that since the 
obstructions are of topological origin the existence of exponentially localized WF is stable against 
perturbations. Indeed along these lines it has been possible to prove the existence of (generalized) 
WF for a variety of nonperiodic systems [2TJ, [55] , [23] , 0) [IS]- Since in the periodic case 
the obstructions to the existence of exponentially localized WF are absent [13] , [14] , [15] in one 
dimension it was naturally to conjecture [H],[24] that in one dimension WF exist for all isolated 
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bands irrespective of periodicity properties. 

The first problem to be solved was to find an alternative definition of WF. The basic idea goes 
back to Kivelson [35] , who proposed to define the generalized WF as eigenfunctions of the "band 
position" operator. To substantiate the idea one has to prove that the band position operator 
is self-adjoint, has discrete spectrum and its eigenfunctions are exponentially localized. For the 
particular case of a periodic one dimensional crystal with one defect Kivelson proved that the 
eigenfunctions of the band position operator are indeed exponentially localized and asked for a 
general proof. In the general case, by a bootstrap argument, Niu [24] argued that the eigenfunctions 
of the band position operator (if they exist) are at least polynomially localized. In full generality 
the fact that for all isolated parts of the spectrum the band position operator is self-adjoint, has 
discrete spectrum and its eigenfunctions are exponentially localized has been proved in [26j . 

In this paper we extend the results in [26] to quasi one-dimensional systems i.e. three dimen- 
sional systems for which the motion extends to infinity only in one direction. In addition we add 
the result (which is new even in the strictly one dimensional case) that (see Theorem [5] below for 
details) the "density" of WF is uniformly bounded. While the main ideas of the proof are the same 
as in [26] there are major differences both at the technical and physical level. In particular for 
quasi one dimensional systems with screw symmetry the constructed WF inherits this symmetry 
a property which is very useful in computational applications. Finally let us point out that as in 
the periodic case, generalized WF defined as eigenfunctions of the band position operator have 
very nice properties e.g. they are (up to uninteresting phases) uniquely defined and for real (i.e. 
time reversal invariant) Hamiltonians they can be chosen to be real functions and this solves for 
the general quasi one dimensional case the "strong conjecture" in Section V. of [3]. As for their 
exponential localization we have the following "optimality" result (see Proposition [3] for a precise 
statement) which seems to be new even in the one dimensional periodic case: the eigenfunctions 
of the band position operator have the same exponential localization as the integral kernel of the 
projection operator of the corresponding band. 



2 The results 

Consider in L 2 (R 3 ) the following Hamiltonian describing a particle subjected to a scalar potential 
V: 

H = P 2 + V, P = -£V, sup / |F(y)| 2 dy<oo (2.1) 

x£l 3 J\x-y\<l 

which, as is well known (see [27]), is essentially self-adjoint on C§°(1R 3 ). We have already said in 
the introduction that we are interested in potentials V which tend to zero as the distance from 
the Ox\ axis tends to infinity. Let us now be more precise. The notation x = (xi,xj_) will be 
used throughout the paper. For any R > 0, define: 

I V (R)~ sup f \V(y)\ 2 dy. (2.2) 

x 1 GR,\x ± \>R J \x-y\<l 

The decay assumption for V will be: 

lim I V {R) = 0. (2.3) 

R — >oo 

It is easy to see that [0, oo) c cr(H) (using a Weyl sequence argument), thus the only region where 
H might have an isolated spectral island is below zero. Now suppose that cr is such an isolated 
part of the spectrum and define: 

- E+ := sup{£ : E G cr } < 0. (2.4) 

If r is a positively oriented contour of finite length enclosing a®, then the spectral subspace 
corresponding to <7o is: 

K := Ran(Po), Pq = ^~ I (H - z)- l dz. (2.5) 
2n J r 
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At a heuristic level, due to the fact that the wave packets from JC cannot propagate in the 
classically forbidden region (see (|2.4p and (|2.3|l ). at negative energies the motion is confined near 
the Ox\ axis, i.e. the system has a quasi one dimensional behavior. 

2.1 The technical results 

The following proposition states the "localization" properties of P . On one hand, this give a 
precise meaning to the previously discussed quasi one dimensional character, and on the other 
hand it provides some key ingredients to the proof of exponential localization of eigenfunctions of 
the band position operator. 

Let a € R, and let (X\^ a ) be the multiplication operator corresponding to: 

5q (x) := V(zi-a) 2 + l, (2.6) 
and (X±) be multiplication operator given by: 

ffx(x) := vIxTFTT. (2.7) 
Proposition 1. There exist an > 0, a>± > 0, M < oo such that: 

sup || e Q ii <x ii'» > P e _Q|l<X|l ' a> ||< M, and (2.8) 

|| e *±(X x )p oe ax(X ± ) ||< M ^ (2.9) 

The proof of Proposition!]] will also give values for ay and aj_. In particular a± can be any number 
strictly smaller than WE+. 

We now can formulate the main technical result of this paper. To emphasize its generality we 
stress that its proof only uses the decay condition (|2.3p and the existence of an isolated part of 
the spectrum satisfying (|2.4[) . 

Theorem 2. Let Xh be the operator of multiplication with x\ in L 2 (R 3 ) and consider in K, the 
operator 

X {[ := PoX l{ P (2.10) 

defined on X>(Xy) = X>(Xy) D K. Then 

i. Xii is self-adjoint on T>{X); 

ii. X\\ has purely discrete spectrum; 

ill. Let g G G := cr(Xii) be an eigenvalue, m g its multiplicity, and {W g j}i<j< m m orthonormal 
basis in the eigenspace of X corresponding to g. Then for all fi € [0,1], there exists Mi < oo 
independent of g, j and (3 such that: 

f e 2(l-/3)«||ki-fll e 2/3a : x|xx|| W9j ( x )|2 dx < Mi) (2.11) 

where ay and a± are the same exponents as those provided by the proof of Proposition^ 

iv. Let a € R and L>1. Denote by N(a,L) the total multiplicity of the spectrum of Xu contained 

in [a — L,a + L\. Then there exists M 2 < oo such that 

N(a,L)<M 2 -L. (2.12) 

Finally, we turn to the question of optimal localization properties of our Wannier functions. 
Theorem [5] provides an optimal exponential decay on the transverse direction, but in the parallel 
direction it only implies a decay which is bound by the maximal decay of the resolvent in the 
gap. The conjecture on optimal exponential decay, as stated in Section V of [3], is whether the 
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W g> j 's have the same exponential decay as the integral kernel Vq (x, y) of Pq (which can be larger 
than the maximal decay of the resolvent in the gap; we are indebted to one of the referees for 
pointing this to us). Concerning this issue, we have the following result showing the optimality of 
the "parallel" decay of W g j at the exponential level. 

Proposition 3. Assume that for all a < ao we are given an a priori bound 

sup || e Q(x i!- > Foe" Q<X|l '" > ||< oo. (2.13) 

Then for all a < ao there exists Mi (a), independent of g and j, such that 

I ^\*i-g\\w gd (x)\*dx<Mi(a). (2.14) 

JR3 



Remark. Here ao is the "exact" exponential decay of Vo(x,y). In certain particular periodic 
cases one might obtain a power-like asymptotic behavior of e a °^ Xl ~ yi 'Vo(^, y) in the variables 
x\,y\. We cannot say anything about an eventual asymptotic behavior of e a °' Xl ~ s 'W gi j(x). But 
due to the generality of the setting, we consider our result to be optimal. 

2.2 Further properties of the Wannier basis 

We come now to the case when V (hence H) has additional symmetries. The point here is that 
although the Wannier functions are not eigenfunctions of H, one would like them to inherit in 
some sense the symmetries of H . The reason is that usually the Wannier basis is used in order to 
write down an effective Hamiltonian in /C, and one would like this effective Hamiltonian to inherit 
as much as possible the symmetries of H . 

First we comment on time reversal invariance. Since V(x) is real, H commutes with the anti- 
unitary operator induced by complex conjugation. It follows (see (|2.5p ) that Pq and Xn are also 

real, thus the eigenfunctions of X\\ can be chosen to be real. Hence Theorem [2] provides us with 
a Wannier basis which is time reversal invariant. 

Second we consider the so called "screw-symmetry" along the 0a;i-axis, of much interest in 
the physics of carbon nanotubes. Namely, writing 

x_L = (r,0), r > 0, 6 € [0,2tt), (2.15) 

one assumes that for some 9q <G [0, 2ir) we have: 

V(x 1 ,r,9) = V(x 1 + l,r,9 + 9 ). (2.16) 

Here + 9o has to be understood modulo 2n. Defining the screw-symmetry operators T®° by: 

(T^f)( Xl ,r,0) :=f( Xl -n,r,e-ne ), (2.17) 

one has a (unitary!) representation of Z in L 2 (M 3 ). Taking into account (|2.16p and the fact that 
[— A,T^°] = (use cylindrical coordinates to prove this), one obtains: 

[HX°}=0, (2-18) 

and then from functional calculus and (|2.5|) : 

[P 0) T*°]=0. (2.19) 

In particular, this implies that the family {T^°} ng 2; induces a unitary representation of Z in JC. 
Moreover, from (|2.10[) and l2.19| ) one obtains: 

[T%°,X {] }=nT%°. (2.20) 
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Let p < oo be the number of eigenvalues of X\\ in the interval [0, 1), and let {gj}^ =1 be the distinct 
eigenvalues (each with multiplicity rrij < oo). We have: 

X\\W gj , ai = gj W gj<aj , a, = l,2,...,m gr (2.21) 
From (|2.20l ) and (I2.2ip one obtains that for all gj, ctj, n e Z: 

X^W g} , aj - ( 9j + n)T°°W gj , aj . (2.22) 

Conversely, for every other g 6 o^Xi), choose an eigenvector W g . We can find n € Z such that 
.g + n G [0, 1). Since X||T^°W fl = (.g + ri)T%°W g , it means that 3 + n must be one of the g/s 
considered above. Therefore we proved the following corollary: 

Corollary 4. The spectrum of Xn consists of a union of p ladders: 

G = lf j=1 Gj, G j = {g: g = gj +n, n e Z}, j e {1, 2, ...,p}, (2.23) 

and an orthonormal basis in K, can be chosen as: 

W n>gj<aj := W gj+n>aj := T°°W gj<aj , (2.24) 
7i e Z, j e {l,2,...,p}, a e {l,2,...,m ff3 .}. 

It is interesting to express the effective Hamiltonian PqHPq as an infinite matrix with the help 
of the Wannier basis. For notational simplicity we relabel the pair (gj, ctj) as I € {1, 2, N c = 
Y^j=i m gj} an d write the Wannier basis as {W n! i}„ eZ . ze{i,2,...,jv c }- Note that N c is nothing that 
the number of Wannier functions per unit cell [0, 1). Let 

h e ^ k {m,n):={W m}l ,HW nik ). (2.25) 

The important fact is that in spite of a rotation with an angle #0 for which it might happen that 
to be irrational, from (|2 . 18[) and (|2.24j) one obtains (with the usual abuse of notation): 

h e l ° k (m,n) = hfa(m-n). (2.26) 

Then a standard computation gives the effective Hamiltonian as an operator in (l 2 ) Nc which is of 
standard translation invariant tight binding type: 

(h e e ° ff ct>)i{m) := h 8 i%{m - n)<p k (n). (2.27) 

k.n 

This is another consequence of the quasi one-dimensional character of the motion for negative 
energies. More precisely, it reflects the fact that for arbitrary values of #0, since is a unitary 
representation of Z, one can still develop a Bloch type analysis but with a more complicated form 
of "Bloch" functions: 

* fc (x) - e lfca;i Ufc (x), Ufc (x) - T°°u k { X ). (2.28) 

However, due to the complicated symmetry of the resulting Bloch functions (which does not allow 
to represent the fiber Hamiltonian as a differential operator on the unit cell with " simple" boundary 
conditions), the analysis gets much harder. The Bloch analysis reduces to the standard one (with 
a larger unit cell) for rational values of 

3 Proofs 

This section is devoted to the proof of Proposition [TJ Theorem [5] and Proposition [31 A certain 
number of unimportant finite positive constants appearing during the proof will be denoted by 
M. 

One of the key ingredients in the proofs is the exponential decay of the integral kernel of the 
resolvent of Schrodinger operators. This is an elementary result in the Combes- Thomas- Agmon 
theory of weighted estimates. We summarize the needed result in: 
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Lemma 5. Let W be a potential such that sup xgR3 /| x _ y |< 1 |W^(y)| 2 ^y < oo. Define K := 
p 2 _|_ J4/( x ) as an operator sum, and let h be a real function satisfying: 

heC°°(m. 3 ), sup{|V/i(x)| + |A/i(x)|} = m<oo. (3.1) 
xe» 3 

Fix z E p{H). Then there exists a z > such that 

\\e a » h {K - z^e-^W < M, (3.2) 
We^PjiK - z)~ 1 e~ a ' h \\ < M, (3.3) 

where Pj = -i^,j€ {1,2,3}. 

Without giving the details of the proof of Lemma El for later use we write down a key identity 
in (13.51): under the condition 



l+a z {±iP- \7h±i\7h-P - a z \Vh\ 2 )(K - z)- 1 invertible (3.4) 

one has 

e ±a * h (K - z)- l e T °" h (3.5) 
= {K- z)- 1 ^ + a z {±iV ■ Vh ± iVh ■ P - a z \Vh\ 2 ){K - z)" 1 ]" 1 . 

Then (|3.4[) holds true if for example a z > is small enough. 
3.1 Proof of Proposition [1] 

Take T in (12.51) a contour of finite length enclosing gq and satisfying 

dist(r, a(H)) = ~ dist(a , a(H) \ a ). (3.6) 

Then since \Vg a \ < 1, |Ag a | 2 < 2, the estimate (12.81) follows directly from Lemma [S]by taking ay 
sufficiently small such that for all z E T: 

\\a l{ (iP • V. 9a + iWg a • P - a,, \Vg a \ 2 )(K - 1 1| < b < 1. 

We now prove (EH). If i? > 0, define: 

= -A + (1 - X r)V, (3.7) 

where 

v M-i 1 for ^ 

Xfl(x) -\ for |xj.|>JZ • (3 - 8) 

From (l2~3l) it follows that 



lim inf a(H R ) = 0. 

R— >oo 

In particular, for sufficiently large i?, (i?R — z) 1 i s analytic inside T. Since H — Hr — xhV> then 
using resolvent identities we obtain: 

(H - z)- 1 = (H R - z)- 1 (3.9) 
- (H R - z)- x X rV{Hr - z)- 1 + (H R - z)- x X rV{H - z^xrV^Hr - z)~ x . 

From ()2.5|) . (|3.9j) and the fact that (_ffr?, — z) _1 is analytic inside T one has 

P ° = h J r {HR ~ z ^ X ^ y{ - H - ^XRViHn - z)-\ (3.10) 



sup f \(e a9± X BV)(y)\ 2 dy < oo. (3.11) 

xeR 3 J|x-y|<l 



Notice that for all a > 0: 



Take now a± > such that 13. 41 holds true for all z € T, K — Hr, h = g± and a z = a±. That is 
let us suppose that 

1 + a±(±iP ■ \7g ± ± iV.g ± ■ P - a ± \X7g ± \ 2 )(H R - z)" 1 is invertible (3.12) 
uniformly on T. Then we can rewrite Pg as: 

P = e - a±{x ^ 

" s «l5l 

-(XX) 



e^XRV(H-z)- 1 ] [e^XBViHR-z)- 1 ] (3.13) 



[l + a±(-iP ■ X7g ± - i\7g ± ■ P - a ± \X7g ± \ 2 )(H R - z)- 1 ] 1 dz\ e~ a± 

Due to (|3. 1 1|) the operator under the integral sign is uniformly bounded in z and the proof of 
Proposition [1] is finished provided we can show why we can choose a± as close to as we 

want. The argument is as follows. Choose < aj_ < J E + . Choose a contour T which is very 
close to co, at a distance 6 > 0, infinitesimally small. Using the spectral theorem (or in this case 
the Plancherel theorem), there exists 6 small enough such that the following estimates hold true: 

sup II (P 2 - z)^ 1 II < const, sup max ||P 7 -(P 2 - z) _1 || < const. (3-14) 

zer »erje{i,2,3} 11 

Hence we can find 8 small enough and R large enough such that the operator in (|3.12ll is invertible 
if 

1 + a ± (±iP ■ X7g ± ± iVg± ■ P - ajV^XP 2 - ^(z)y 1 is invertible (3.15) 
uniformly on T. Now the operator in (|3.15p is invertible if 

1 ± ia ± (P 2 - »(z))~'(P • Vfe + V/i ■ P)(P 2 - 
- ai(P 2 - 3f?(z))-5| V /i| 2 (P 2 - SR(z))-2 (3.16) 

is invertible (by a resummation of the Neumann series and analytic continuation). Now assume 
that uniformly on T we have: 

< «i(P 2 - 5i(z))^|V/*| 2 (P 2 - < -|L < 1, 

which can be achieved if a\ < E + and 5 is chosen to be small enough. Define 

S:=(l- a 2 ± (P 2 - 3f?(z))^|V/i| 2 (P 2 - K(.z))-*) _7 , 

and 

T = T* := S(P 2 - 3?(z))"5(P • V/i + V/i • P)(P 2 - 
Then the operator in (|3.16l) is invertible if 1 ± ia±_T is invertible, which is always the case: 

(1 ± iaxT)- 1 = (1 =F ia±T)(l + a^T 2 )' 1 . 

Therefore Proposition [1] is proved. □ 
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3.2 Proof of Theorem 

Proof of (i). First we recall an older result (see e.g. [251 HI US] ) > according to which the commutator 
[Xi,Pq] defined on T>(X\\) has a bounded closure on L 2 (R 3 ). We seek an approximate resolvent 
of Am by defining for [i > the operator 

R±„. = Pv(X\\ ± in)- 1 P - (3.17) 

Since one can rewrite R± /1 as 

i? ±M = (X || ± ^M)- 1 Po + ± i^y 1 ^ , P ] (X|, ± i/^Po 
it follows that R± f _ l IC c -D(Xii) and by a straightforward computation (as operators in IC) 

(X,| ± z M )P± M = P (X|| ± i/x)P (X|| ± iM) _1 Po = Ijc + (3-18) 

with 

A ±fl = P [X||,P ](X|| ±i t i)- 1 P . (3.19) 
Since [Xy,Po] is bounded and ||(Aj| ± i/x)" 1 1 < —, it follows that for sufficiently large fi: 

U±,\\ < \. (3.20) 

Then again as operators in IC: 

(X±i f i)R ±ti (l K + A ±tl )- 1 = l K (3.21) 

This implies that X ± i/i is surjective on R± ij ,{1k + A±^)~ 1 JC C D(X) . By the fundamental 
criterion of self-adjointness [27] X is self-adjoint in K. on D(A). In addition, from (|3 . 2 1 [) one 
obtains the following formula for the resolvent of X\\ : 

± z/i)- 1 = R ±I1 {1 K + A ±ll )-\ (3.22) 



Proof of (ii). We will show that R± ll is compact in K, which implies (see p. 22^ 1 that Xt\ has 
compact resolvent, thus purely discrete spectrum. Consider a cut-off function 4>n which equals 1 
if |x| < TV and is zero if |x| > 2N. For N > 1 we can decompose: 

R ±tl = P (X|| ± i/x)~VivPo + P (X|| ± - MPo- (3.23) 

Writing 

^ N P = {0iv(P 2 + 1) _1 }{(P 2 + l)Po} 

we see that 4>nPo is compact (even Hilbert-Schmidt) in L 2 (R 3 ) (the first factor is Hilbert-Schmidt 
while the second one is bounded). Now if < a is small enough, we know that e ag± Po is bounded 
(see ([2310 . Since 



lim \\(Xu ±ifi)- 1 (1 - </) N )e- agx II = 0, 



we have shown: 



lim 

N— >oo 



P± M -Po(X|i ±iM)"VwPo 



0. 



thus P± M equals the norm limit of a sequence of compact operators, therefore it is compact. Ac- 
cordingly, since the self-adjoint operator Xi\ has compact resolvent it has purely discrete spectrum 

= o-disc{X\\) =: G, (3.24) 
and the proof of the second part of Theorem [2] is finished. 
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Proof of (iii). Now we will consider the exponential localization of eigenfunctions of X\\. Let 
g £ G be an eigenvalue, m g its multiplicity and W g j, 1 < j < m g be an orthonormal basis in the 
eigenspace of X\\ corresponding to g. We shall prove that uniformly in g and j 

\\e a " {x ^s } W g ,j\\ < M and (3.25) 
\\e a±{x±) W gj \\ < M. (3.26) 

Taking (|3.25[) and (|3.26p as given, one can easily obtain (|2.11j) by a simple convexity argument: 
the function f(x) = a 1 ^ x b x ; a, b > is convex on M., and for < [3 < 1 one has: 

e 2a|| S a(x) + ^ _ fl e 2*_Lgx > e 2(l-/3) Q | |ga (x) e 2f3 a±g± ^ ^.21) 

which together with (|3~2"5]) and (pT2l)|) it proves ([2"TTj) with M x = M 2 . Since (|3T2l)]) follows directly 
from (|23|1 and W g>j = PoW gJ we are left with the proof of (pHS)) . 

Although the proof of (|3.25p mimics closely the proof in the one dimensional case [55] , we give 
it here for completeness. In order to emphasize the main idea of the proof let us remind one of the 
simplest proofs of the exponential decay of eigenfunctions of Schrodinger operators corresponding 
to discrete eigenvalues (assuming that the potential V is bounded and has compact support). 
Namely assume that for some E > we have (—A + V + E)fy = 0, which can be rewritten as 

* = -(- A + E)- 1 Vy. (3.28) 

Since for \a\ < y/E, e Q N (-A + E^e' 01 ^ and e a ^V are bounded: 

* = -e"" 1 ' 1 {e l - l(-A + ^)- 1 e -0 l-lJ(e l , IV r )* 

which proves the exponential localization of \&. The main idea in proving (|3.25p is to rewrite the 
eigenvalue equation for Xn in a form similar to (|3.28p and and then to use (|2.8p . 
Let us start with some notation. If b > (sufficiently large) and a E K, define: 

/a,6(x):=&/(^p) (3.29) 

where / is a real Co°(R) cut-off function satisfying < f(y) < 1 and 

1 for | v | < i 



for 1 2/| > 1 



f(y) = 

Define the function h a ,b by: 

h atb (x) :=xi-a + ifaA x )- {3.30) 
Note that by construction, h aj b only depends on xi, and obeys: 

IVfcWI > \- (3.31) 
Moreover, its first two derivatives are uniformly bounded: 

sup supsup{|V/i aj b(x)| + |A/i ajb (x)|} = K < oo. (3.32) 

xSR 3 aGRb>l 

The eigenvalue equation for W g j reads as -Po(^|| — 9)PoW g .j — 0. Using (|3.30p it can be rewritten 
as: 

P h g>b P W g!j = iP f gib P W gtj . (3.33) 

We now prove that Pohg^Pa is invertible. Like in the proof self-adjointness of Xn we compute 

P o h- t lP o P o h g , b P = 1 K + Pohg l [Po, h g ,b] Po- (3.34) 
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The key remark is that [Pq, is bounded. Indeed we have the identity: 

[Pa, h g . b ] = --!- f(H- z)- 1 {P • Vh g , h + Vh g . b -P}(H- z)- x dz 

= f(H- z)- 1 {-iAhg b + 2X7h g b -P}(H- z)- x dz. (3.35) 

It follows that [Pq, h g: b] is uniformly bounded in g G K and b > 1 (see (|3.32p ). Taking into account 
(|3.3ip one obtains that the operator 



B g<b =P h-l[P ,h gtb ]P : K^K 



satisfies 



v 9 A\ < o 



(3.36) 



(3.37) 



if b > bo for some large enough b < oo. It follows that l + B gb is invertible and then the eigenvalue 
equation (see (|3.33p , (|3.34p and (|3.36p ) takes the form 



g^-' UJ g,t>* U" g,j (3.38) 

which is the analog of ()3.28j) . By construction (see the definition of f g .b in (|3.29|) 1: 

||e a H <X| l' 9> / 9 ,6|| < 6e a » (b+1) . 

Moreover, 

e a \\( x \\. g ) P h' g lP e~ a ^ x ^^ = |e Q ^ x ii' 9 ^ P e~ a ^ x ^^\ h~\ ie a ^ x hg> P e~ a ^ x ^^\ 

is bounded due to (|2.8|1 . Thus the only thing it remains to be proved is the existence of a 6 large 
enough such that the following bound holds: 



sup 



< oo. 



(3.39) 



Using the Neumann series for I 1 + B g b ) , it follows that it suffices to prove that 



lim sup 



; a «( X \l,g) B , e - a \l( X H,g) 



(3.40) 



Since (see (|3.3ip ) lim^oo \\h g l\\ — (uniformly in g 6 M), for (|3.40p to holds true it is sufficient 
to show: 



sup 



;°U<*II,S> [P ,h gb ] e- Q ll (X ll^) 



< const. 



(3.41) 



But this easily follows from (|3.35|) , (|3.32p . (|3.2p and (|3.3[) where we take K = H, a z = an and 
h = g g . The proof of (iii) is concluded. 

Proof of (iv). We start with a technical result: 

Lemma 6. Fix < ctj_ < E + . Then there exists a bounded operator D such that 

P = e - a ^ x ^ (P 2 + l^D (3.42) 
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Proof. We use the notation and ideas of Proposition [1] and we rewrite Pq in a convenient form. 
First, for R > we have 

(H - z)- 1 = (H R - z)- 1 - (H R - z)- l XR V{H - z)-\ 

Second, choose T close enough to ao an d R large enough, such that (Hr — z)^ 1 becomes analytic 
inside T and (|3.12[) holds true for all z € V. Then we can write: 

Po = -e- a ^ x ^ 'JL /" (3.43) 

- z)- 1 !! + a±(iP ' V.g± + iVs± • P - a ± \Wg±\ 2 )(H R - z)- 1 r 1 e a ^ XR v(H - z^dz. 

Now by the closed graph theorem we have that (P 2 + 1)(H R + is bounded (here R is large 
enough such that (— oo, —1/2) C p(H R )), and together with the spectral theorem: 

sup\\{P 2 + l)(H R -z)- 1 \\ <oo. 
zer 

Use this in (|3.43[) and we are done. □ 

We now have all the necessary ingredients for proving the last statement of our theorem. For 
every L > and agl, denote by XL,a the characteristic function of the slab {x : \x\ — a\ < L}. 
Then define the operator B := XL.aPo- Using (|3.42|) let us show that B is Hilbert-Schmidt, and 
moreover, uniformly in a G M we have: 

\\B\\l<M-L, (3.44) 

for some M < oo. Indeed, since B = XL,a,e~ a±i ' X± ' > [— A + 1)~ 1 D, a direct computation using the 
explicit formula for the integral kernel of the free Laplacian gives: 

\\XL, a e- a±{Xx) (P 2 + l)" 1 ^ < const • L, 

It follows that the operator XL,aPoXL,a = BB* is trace class and 

\Tr(XL, a PoXL,a)\ < \\B\\l <M-L (3.45) 

for some M < oo independent of L and a. 

Now let Pg' a be the orthogonal projection onto the subspace spanned by those W g j for which 
g E [a — L,a + L\: 

rn g 

P o' a -= E E<-.^j)W w -. (3.46) 

\g-a\<Lj=l 

We can choose A sufficiently large such that (|3.25p implies: 

|IU 9J (x)| 2 dx<i (3.47) 

\x x -a\>A z 

uniformly in a and g G [a — L, a + L\. Since P > P^ a , from (I3~l5l) one obtains: 

M ■ (L + A) > Tr( X L+A,aPoXL+A:a) > Tr( X L+A,aPo' a XL+A:a) 



E E / XL+A,a(x)|IU SJ (x)| 2 dx 



^ E fl\ = \ N ^ L )i ( 3 - 48 ) 

\g-a\<Lj=l 

where in the last inequality we used (|3.47p . In particular, if L > 1, then uniformly in a G R we 
have 

N(a,L) < 2M ■ (1 + A)L 

and the proof is finished. 
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3.3 Proof of Proposition [3] 

The only thing we have to prove is that (|3.41[) holds true for ay replaced by any a < ao, where 
«o is the a-priori given, "exact" exponential localization. 

We introduce the multiplication operator given by {e a '' _t '/}(x) := e a ' :El ~ t '/(x). We start by 
noticing that due to the bound e ±a W s2 + 1 -\ s \) < e Q we can replace (|2.13[) with: 

sup ||e a l-*li%e- a l-*l || < oo. (3.49) 

ten 

The same replacement can be done in (13.41[) . Now the integral kernel A(k, y) of the operator 
A := e a \- g \ [P ,h g>b ] e'^-^ equals 

A(x,y) =V (x,y)e a ^-^-^(h 9 Ay) ~ h g A*))- (3-50) 

We consider A as an operator on L 2 (M 3 ) = (J) peZ L 2 ([p,p+ 1] x K 2 ). Let \ p be the characteristic 
function of the slab [p,p+ 1] xl 2 . We have that A pp * :— XpA-Xp 1 1S a bounded operator between 
L 2 ([p' \p' + 1] x M 2 ) and L 2 (\p,p + 1] x M 2 ), and we can write A — {A pp /} PiP / eZ . We will bound 
the norm of A with a Schur-Holmgren type estimate (see below Lemma [7]) : 




HAH < sup^H^II Up£|LV|l ■ (3-51) 
Y^pez J \^p>& ) 

For < Xi,yi < 1, the kernel of A pp i can be written as: 

-V( x >y) = +p,*±;vi + v' ,y±)e^ xi+p - a ^ yi+p '-^(h g>b ( yi + P ') - h g>b ( Xl + P )) 

= P (xi + P ,x±;yi +p',yi_)e a ^+P-a\-^+p'-3\)^ b (p') - h g , b {p)) 

+ V (x 1 +p,x ±;yi + P ',y±)e a ^+P-^-^ + P'-^(h g , b ( yi + p') - h g , b (p')) 

+ V Q { Xl +p, X± ; yi + p', yi_)e a ^+P-a\-\y^p'-9\){-h a ^ Xl + p) + h g , b (p)) 

=: A%, (x, y) + A$ (x, y) + A p 3 p \ (x, y). (3.52) 

The last two kernels can be analyzed with the same methods as the first one, thus we only estimate 
the norm of A pp \ . The crucial observation is that we can write this operator as a product of three 
operators having the corresponding kernels: 

A pp ](x,y) = e «(l*i4*- fl |-|p- fl |) 

. e «Up-g\-tf-g\) Vo ( Xl +p^ riVl +p',y x )(h g , b (p') ~ h g , b (j>)) 

. e -<*(\vi+p'-g\-\p'-g\) _ (3.53) 

The kernel in the middle corresponds to the operator XpPoXp times some coefficients depending 
on p,p' . 

Using the triangle inequality to bound the exponentials, and (|3.32p in order to write \h gtb (y) — 
h g A x )\ ^ K\ x i - ViV we nave: 

H^SlI < Ke 2a e a \P-P'\\p~p>\ ■ \\XpP x'p\\- 
Using t — p' and (a + ao)/2 in (|3.49| we obtain 

thus 

\\A pp \\\ < C'\p - p'\e-( a °- a ^~r'\/ 2 
which is summable in the sense of (|3.51|) . The same strategy can be applied in the case of A pp , 

(3) 

and A pp , . The last thing to be done is to prove the Schur-Holmgren estimate: 



13 



Lemma 7. The estimate (|3.5ip holds true. 

Proof. Let ip € L 2 (M 3 ) with compact support and \ \ip\ \ = 1. We write: 



PV>l| 2 =£ll^H 2 . (3.54) 

pez 



But 



IIxp^II < £ Vll^'H Vll^VH Hxp^H ^ i E II^VII \ \ E Ppp'H 11x^1 

p'GZ [p'eZ J [p'GZ 

^{^pEii^h} |eii^'IIII^'^II 2 | ( 3 - 55 ) 

where in the second inequality we used Cauchy-Schwarz with respect to p' . Introduce this in (|3.54[) 
and the bound follows after the use of J2 P 'ez IIXp'V'll 2 = 1- d 
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